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Morphology of the complex HI gas distribution can be quantified by statistics like the Minkowski 
functionals, and can provide a way to statistically study the large scale structure in the HI maps 
both at low redshifts, and during the epoch of reionization (EoR). At low redshifts, the 21cm emis¬ 
sion traces the underlying matter distribution. Topology of the HI gas distribution, as measured by 
the genus, could be used as a “standard ruler”. This enables the determination of distance-redshift 
relation and also the discrimination of various models of dark energy and of modified gravity. 
The topological analysis is also sensitive to certain primordial non-Gaussian features. Compared 
with two-point statistics, the topological statistics are more robust against the nonlinear gravita¬ 
tional evolution, bias, and redshift-space distortion. The HI intensity map observation naturally 
avoids the sparse sampling distortion, which is an important systematic in optical galaxy survey. 
The large cosmic volume accessible to SKA would provide unprecedented accuracy using such 
a measurement. During the EoR, topology can be a powerful and intuitive tool to distinguish 
among the different evolutionary stages of reionization, where the ionized regions make up a sig¬ 
nificant fraction of the volume. Eurthermore, it can also discriminate among various reionization 
models. The genus curves evolve during cosmological reionization, and for different reionization 
scenarios, the topology of the HI gas distribution can be significantly different even if the global 
ionization fractions are the same. It can provide clear and intuitive diagnostics for how the reion¬ 
ization takes place, and indirectly probes the properties of radiation-sources. In this brief chapter 
we will describe the scientific background of the topology study, and forecast the potential of the 
SKA for measuring cosmological parameters and constraining structure formation mechanism 
through the study of topology of HI gas distribution. 
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1. Topology study of HI distribution 

Topology has been introduced in cosmology to describe the properties of the large-scale struc¬ 
ture (LSS), and to test the non-Gaussianity of the primordial density fluctuations (Gott et al. 1986; 
Hamilton et al. 1986; Park et al. 2013). Application of LSS topology study has been extended 
to measurement of cosmological parameters and constraining galaxy formation mechanism (Park 
et al. 2005; Choi et al. 2010, 2013). The LSS topology is also a cosmological invariant that can be 
used to reconstruct the expansion history of the universe (Park & Kim 2010; Zunckel et al. 2011; 
Wang et al. 2012; Speare et al. 2013; Blake et al. 2014). Mathematically, geometry of the excursion 
regions with density above a threshold can be characterized by the Minkowski Functionals (MFs) 
' (Mecke et al. 1994; Pratten & Munshi 2012; Hikage et al. 2006, 2008). The genus, related to 
one of the Minkowski functionals, is given by the number of holes minus the number of isolated 
regions in the iso-density contour surfaces of the smooth density field, and is a measure of topology 
that quantifies the connectivity of the contour surfaces (Park et al. 2013). All the MFs are analyt¬ 
ically known for Gaussian fields, the genus per unit volume as a function of the threshold density 
V = 5/a is given by 


G{v)=A{l-v^)e-^"l\ A = (1.1) 

The topology of the isodensity contours is insensitive to the systematic effects such as redshift- 
space distortion, non-linear evolution, sparse sampling (Park et al. 2005; Park & Kim 2010). 
Therefore, any deviation from the Gaussian prediction is evidence for non-Gaussianity, and can 
constrain the mechanism for primordial non-Gaussianity. The deviation can be caused not just by 
the three-point function but by all the high-order moments of the density field. Therefore, the study 
of non-Gaussianity using MFs is complementary to the approach of using three-point correlation 
function or bispectum, and is in principle able to detect general forms of non-Gaussianity. The 
Minkowski functionals measured from the HI 21cm map can characterize the complex distribution 
of the neutral hydrogen, both in the large-scale structure at low redshifts, and during the epoch of 
reionization (EoR) at high redshifts. 

1.1 Primordial non-Gaussianity 

The primordial perturbation may deviate from the Gaussian random field during the cosmo¬ 
logical phase transition triggered by the spontaneous breakdown of symmetry creating topological 
defects such as the cosmic monopole, string, wall, and textures. Various non-standard inflationary 
scenarios which violate the single-field and slow-rolling condition may also produce significant 
primordial non-Gaussianity. It is usual to parameterize the non-Gaussian features generated by the 
inflation model as <P(x) = -|-/nl “ (^g)] where (pQ is the Gaussian distribution of 

potential. For this form of primordial non-Gaussianity, the genus curve deviates from the Gaussian 

*In three dimensions, the four Minkowski functionals for the isodensity contours are the volume fraction Vo(v), 
surface area H(v), mean curvature V 2 (v) and Euler characteristics V 3 (v) = ;t(v), which is related to the genus by 

;^ = 2”2g. 
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model, and up to second order perturbations it is given by Hikage et al. (2006, 2008). 


A(v) = 5 


G 


= -e X /nl 


(s^ - s'^mvf )+- s;s>)H,(v) 


ao 


( 1 . 2 ) 


where is the skewness parameter (Matsubara 2003). The measured genus curve can then be 
used to constrain //^l- However, even if the non-Gaussianity can not be parameterized in this form, 
it would still affect the genus curve, so the topology can be used to probe more general forms of 
primordial non-Gaussianity. 


1.2 Dark Energy and Modified Gravity 

The LSS topology is relatively stable with the growthing of structure. Eq. ( jl.lD shows that A 
measures the slope of the power spectrum around the smoothing scale Hq, and for linear growth A 
would be conserved (Park et al. 2005), so as long as the same smoothing scale is used, the same 
comoving volume at different redshifts would enclose nearly the same amount of structures, as 
illustrated in the left panel of Fig.|I] This fact can be used to measure the redshift-distance relation 
r(z) (Park & Kim 2010). Since the topology of the structure is not scale-free, the genus enclosed in 
a wrongly sized volume and smoothed with a wrong scale would lead to deviation from the actual 
one. 


M(z,Rg,y)R}}.y =^x{Rg,x)Rg,x^ (1-3) 

where R is the volume, Y is the adopted cosmology parameters while X is the true cosmology. The 
smoothing scale Rq for different cosmologies are related by {Rg,x/Rg,y)^ = {Da/^)x/{Da/^)y: 
where Da,H are the angular diameter distance and Hubble expansion rate respectively. Utilizing 
this effect, topology of the large scale structure can serve as a standard ruler in cosmology. This 
technique has been applied to constrain dark energy model (Blake et al. 2014; Zunckel et al. 2011). 

For modified gravity models, the growth rate is different from the general relativity, and in 
some cases it is scale-dependent. This will induce changes in the genus curve. In the right panel of 
Fig.|]the redshift evolution of the genus amplitude for the f{R) theory and some phenomenological 
modified gravify models are shown. Here fhe modified gravify model paramefer Bq is fhe presenf 
day value of fhe funcfion B{a), which is fhe square of fhe Compton scale given by B{a) = 
where /r and /rr are fhe firs! and second derivafives of f{R), wifh R being fhe Ricci scalar. 

So fhe topological measuremenfs can also be used to consfrain fhe modified gravify models 
(Wang el al. 2012). Compared wifh direcf measuremenf using fhe growfh faclor, fhe fopological 
measuremenf may be more robusl, as fhey suffer less from fhe effecf of bias, redshiff distortion, 
and non-linearify. 

1.3 Reionization 

The topology analysis provides a direcf and sensitive probe for fhe delailed process of cos¬ 
mic reionizafion (Fee el al. 2008; Hong el al. 2014). Using fhe genus of HI densily conlours as 
a quanfilalive measure of fopology, fhe reionizafion process of fhe inlergalaclic medium (IGM) 
can be divided into four dislincl topological phases for fhe slandard scenario of reionizafion: (1) 
pre-reionizalion, before fhe formation of any HII bubbles, fhe HI disfribufion reflecls fhe primodial 
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Figure 1: Left; The number of structures as cosmic standard ruler. Right; the redshift evolution of genus 
amplitude for f{R) and phenomenological modified gravity models with different Bq values. 


density fluctuations, the genus curve is consistent with a Gaussian density distribution, and topol¬ 
ogy remains unchanged as the HI density evolves linearly; (2) pre-overlap, characterized by a 
topology dominated by isolated HII bubbles. At the earlier stages of this phase, the number of 
HII bubbles increases gradually with time, resulting in an increase in the amplitude of the genus 
curve, while at the later stages when HII bubbles start to merge, the genus curve goes through a 
turnover and subsequently decreases in amplitude. (3) overlap, during which the amplitude of the 
genus curve drops significantly as the HII bubbles rapidly merge; (4) postoverlap, when the IGM 
is almost ionized, and the evolution of the genus curve is consistent with an decreasing number of 
isolated neutral islands. The measurement of the genus curve can be used to characterize the evo¬ 
lutionary stages of the neutral topology and to identify the redshifts at which the different stages 
of reionization occur, thus distinguishing different reionization scenarios. This in turn provides 
information on the nature and properties of the early luminous sources. 


2. The SKA performance 

The observation of HI emission is presently limited to z < 0.2 by the sensitivity of existing 
telescopes. The SKA, with its great increase in sensitivity, will allow us to observe the HI at much 
higher redshifts, and probe the topology of HI distribution with unprecedented capability. The 
SKA will be built in two stages (Dewdney et al. 2013; Braun), the SKA-1 will consist three arrays; 
SKA-mid, SKA-sur and SKA-low, the features are listed in Table.|T| The SKA2 is to be designed 
in the future, but should be about 10 times larger in mid-frequency and 4 times larger in the low 
frequencies. 

The SKA-low is designed for observation of the high redshifts of the EoR, while SKA-mid for 
the low-to-mid-redshifts of post-reionization Universe, and SKA-sur with its larger FoV is suitable 
for survey of large area. All of these can be used for HI topology study. The topological analysis 
can be applied on data from all HI surveys. We expect that large HI surveys will be one of the key 
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Table 1: The three SKAl arrays. 


Array 

configuration 

relevanf band 

SKA-mid 

190 X 15m dish -i- 64 MeerKAT 13.5m dish 

band 1(0.35-1.05GHz), band 2(0.95-1.76GHz) 

SKA-sur 

60 X 15m dish -i- 36 ASKAP 12m dish (PAF) 

band 1 (0.35-0.5) and band 2 0.65-1.67 GHz. 

SKA-low 

250,000 log periodic anfenna in 866 slafions 

50-350 MHz 


science projects of SKA, being carried out with multiple applications, and the HI topology can be 
studied in conjunction with the other projects. 

A classic approach is the HI galaxy survey where the individual galaxies are observed as 
objects whose HI content exceeds detection threshold: 


Slim — 




sys 


= Nth- 


SEFD 


Ag/yv/Avt v/Avt 


( 2 . 1 ) 


where Nth is a preset threshold value, e.g. 5 or 10. The SKA is potentially capable of surveying 
a billion of galaxies, with SKA-1 it is capable of surveying 10^^^ galaxies(Abdalla & Rawlings 
2005; Myers et al. 2009) at z up to ~ 1.5. But even smaller surveys, e.g. those conducted with 
SKA-1 early science, or for small deep field, would for the first time reveal the HI distribution at 
redshifts much beyond our current limit, and the topological analysis will help us to understand 
quantitatively the characteristics of HI distribution and redshift evolution, as has been down for 
low redshift optical galaxies (Choi et al. 2010, 2013). 

Another approach is intensity mapping (IM): the sky is mapped at low angular resolution, 
such that individual galaxies can not be distinguished, but data of HI large scale distribution can 
still be obtained (Chang et al. 2008). For the EoR experiments, including that of the SKA-low, 
due to the low angular resolution, this will also be the mode of observation. A few dedicated 
experiments, such as the CHIME^ and Tianlai ^ are designed to conduct IM observations. While 
the SKA-mid or SKA-sur are not designed for this, it has been proposed that IM surveys may be 
conducted with the dishes used individually with output autocorrelation, while being calibrated 
using the interferometry (Santos 2014). Whether this would work still needs to be tested in the 
field, buf if successful, if would allow higher sensifivify on inferesfed scales af fhe price of lower 
angular resolufion. 

In Table^ we lisf bofh fhe HI galaxy and IM survey paramefers we adopfed for our forecasl 
on dark energy and modified gravify. For SKA-low, as fhe design is sfill very uncerfain, we do nol 
make a full forecasl here, buf only discuss fhe potential of fhe measuremenf. 


3. Large Scale Structures 

Afler EoR, wilh fhe IGM highly ionized, fhe neufral hydrogen resides mainly in galaxies, 
and fhe HI topology fraces ouf fhe ESS. We can use fhis observation to consfrain fhe primordial 
non-Gaussianify, dark energy, and modified gravity. 


^http://chime.phas .ubc .caj 
^http://tianlai.bao.ac.cn/ 


5 














HI topology 


Xuelei Chen 


Table 2: SKA survey parameters utilized in this paper. 


galaxies surveys 

survey area (deg^) 

integration time (h) 

flux limit (mJy) 

SKA 1-Mid 

1,000(5,000) 

5,000(25,000) 

0.13 

SKA2-Mid 

1,000(5,000) 

5,000(25,000) 

0.05 

SKAl-Sur 

30,000 

7,500(50,000) 

0.4 (0.15) 

SKA2-Sur 

30,000 

7,500 

0.05 

HI intensity mapping 

survey area (deg^) 

integration time (h) 

frequency band (MHz) 

SKAl-Mid-SD 

30,000 

10,000 

350- 1050 

SKAl-Mid-Int 

30,000 

10,000 

350-800 

SKA2-Mid-Int 

30,000 

10,000 

350-950 



Figure 2: The measurement error of /nl from genus curves. We mark the survey volumes of previous 
(SDSS-DR7, BOSS) and forth-coming (DESI & SKA) surveys. 


For primordial non-Gaussianity, Fig.§ shows the /nl constraints from measurement of genus 
curve in the Horizon Run III simulation, which has a volume of 10.8^/j^^Gpc^. We employed a 
Gaussian smoothing of Rq = 22/i^*Mpc to build a density field from subcubes, and compare the 
result for four survey volumes: SDSS, BOSS, DBS and SKA. If the SKA-2 will observe galaxy 
distribution out to z = 3 with a solid angle of 20,000 deg^, then using the genus-curve we can 
measure /nl with an error, cj(/nl) = 20. We should note that this is comparable to the CMB 
limits, and the genus measurement can detect more general forms of non-Gaussianity. 

Using the genus curve as standard ruler, dark energy equation of state parameters wo,Wa can be 
constrained. As shown in Table (^, we consider a relatively deep survey using SKA-Mid compared 
with shallow but wider surveys for SKA-Sur. 

The HI galaxies survey constraints are plotted in Fig. (||), and IM survey constraints are plotted 
in the right panel of Fig.^ With a conservative assumption of the survey parameters that require 
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Figure 3: Constraint on dark energy with HI galaxy survey. Left: HI galaxy survey with SKAl-mid, Right: 
HI galaxy survey with SKA 1-survey. The dashed contours indicate the constraints from Euclid. 


SKA—mid 21cm intensity mapping 




Figure 4: SKA-mid HI intensity mapping. Left: Various comoving scales as a function of redshift. The 
upper red region denotes the quasi-linear regime where the systematics of genus statistics are well-controlled. 
Blue bands are most sensitive scales that would be observed by the core of SKA-mid interferometer, where 
we assume the baseline range from 100m (50m for lighter blue band) to 1000m. The grey region at bottom- 
right corner are scales available with the SKA spiral arm baselines. The black dashed line indicates the 
physical angular resolution of single dish. The solid black line is the effective resolution [R^^^RfreqY^^, 
where we utilize the minimum quasilinear scale, i.e. the boundary of red region, as Rfreq- Right: Constraints 
on dark energy equation of state wq and Wa- In the calculation, ‘varying R’ means the smoothing length is 
selected using smallest quasi-linear scale. The dashed contours indicate the constraints from Euclid. 
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5,000 hours integration time, neither constraint is very stringent for SKAl, but it is independent of 
other techniques, and strong constraints can be achieved with SKA2. As we already know that BAO 
gives better constraints than genus for dark energy, we just give the constraints of our calculation, 
and refer the interested reader to the chapter written by Bull et ah. 

For the IM, there is certain complications. Unlike the BAO technique, where the physical scale 
is fixed, the topological standard ruler measure the scale around the smoothing length. To achieve 
better statistics, smaller smoothing length is desirable, though the systematics such as nonlinear 
growth put limits on the usable scale. In the left panel of Fig.^, we illustrate various physical 
scales for SKA-mid. The quasi-linear region where systematics is well controlled is indicated by 
the upper red region. 

For the single dish observing mode, the angular resolution by the 15m dish is poor (black 
dashed-line), varying from 30Mpc//j at z = 0.5 to above 200Mpc//j at z = 3. However, in the 
line-of-sight direction the resolution is good, the effective 

/?eff=(/?Lg^freq)^/^ (3-1) 

becomes much smaller. Assuming a reasonable 7?freq as indicated by the lower edge of the red 
region, the effective smoothing length varies from 10Mpc//j to 60Mpc//j. The cosmological con¬ 
straints of this observation mode is shown as the blue contour in the right panel. 

If sufficient number of short baselines is available, it may be possible to make an IM survey 
in the interferometric mode. We assume a compact core in the center of SKAmid, the blue belt in 
the left panel illustrates the most sensitive physical scale of this core, where we assume baselines 
distributed from 100m (50m for lighter blue region) to 1000m. With the systematic limit from 
nonlinear scale, one would be able to observe higher redshift volume from z = 0.78 (z = 0.5) at 
much smaller smoothing length and therefore a better statistical accuracy. As shown in the right 
panel, such observation mode improve the constraints significantly. 

For modified gravify models, we consider specifically the f{R) theory with a single model 
parameter Bq characterizing the present value of Compton scale. HI galaxy survey by SKA2 will 
be able to constrain Bq to 5 x 10^^ at 1-a level for a fiducial value of Bq = 10^^, while fhe best 
constraints from intensity mapping is 8 x 10^^. This is much better than the current strongest 
constraints of 1.1 x 10^^ at 95% C.L.. 

4. Reionization 

The genus curve can be used to distinguish the different phases of reionization (Lee et al. 
2008; Hong et al. 2014) discussed in § [L3| . In Fig. || we show the genus curve from an EoR simu¬ 
lation (Trac et al. 2008), for four redshifts during the EoR. As expected, at the highest redshift(z = 
19.7,X, « 0), the genus curve has the shape for that of a Gaussian distribution. The amplitude 
of the curve increases in the pre-overlap phase (z = 9,x, = 0.14), but retains the shape for Gaus¬ 
sian distribution. The shape of the genus curve is drastically different during the overlap phase 
(z = 6.97,X, = 0.65). Einally, the curve returns to the shape of Gaussian distribution in the post¬ 
reionization phase (z = 5.99,x, = 0.99). Erom Eq jl.l| we see the amplitude of genus curve A oc kl, 
ke first increases as the HII bubbles grows in number, and then decreases as the bubbles merge and 
overlap. Measuring k^ from A can provide information on the bubble sizes. 
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V V 

Figure 5: Genus curve and statistical error of 57], at different redshifts from a simulation (Trac et al. 2008), 
the simulation box has a side length of lOOMpc h^*, and smoothing length is IMpc h^* 

These curves are obtained with a smoothing length of lMpc//i. Along with the genus curves, 
the statstical error obtained from a 100Mpc//i simulation box is plotted. The design of the SKA- 
low (Dewdney et al. 2013; Braun 2013) is still very uncertain, but it is agreed that there should be 
a compact core within 1km which provides high sensitivity observation on scales of a few arcmin 
for redshift 10. This is roughly sufficient for the measurement. The observation of SKA-low for a 
deep field should be able fo yield dafa which can be compared wifh such simulations. 

5. Conclusion 

Topology can be very useful at characterizing random fields. The SKA, wifh ifs high sensifiv- 
ify, can observe fhe HI disfribufion much beyond fhe currenf limifs. We can apply fhe fopological 
analysis mefhod fo fhe SKA HI survey dafa. Such analysis can also be used fo fesf or consfrain 
primordial non-Gaussianify, dark energy and modified gravify. The consfrainfs for fhe model pa- 
ramefers from fhe topological mefhod are nof as fighf as fhe fradifional mefhods, e.g., dark energy 
consfrainfs from BAG measuremenfs (See Bull ef al. chapfer) and primordial non-Gaussianify 
from scale-dependenf bias and bispecfrum measuremenfs (See Camera ef al. chapfer), buf if is very 
robusf againsf non-linearify, bias and redshiff disforfion in fhe evolution of large scale sfrucfure. 
During fhe EoR, if can also distinguish differenl evolufion phases and characters of reionizafion 
model. We made forecasls for some of fhese applicafions, in bofh HI galaxy survey and HI in¬ 
tensify mapping survey, and obfained consfrainfs for differenf survey parameters. For SKAl, fhe 
precision of fhe measuremenf wifh fhe genus curve is moderafe, buf we anticipate fhaf fhe addifion 
of ofher Minkowski funcfionals would help improve fhese consfrainfs, and for SKA2 fhe resulf is 
compefifive. More imporfanfly, fhe genus curve is capable of probing non-Gaussian feafures which 
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are not parameterized in the standard form, so with SKA it would be an important tool for probing 
and diseovering non-Gaussianity in unexplored regimes. 
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